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We describe the large-scale collective behavior of solutions of polar biofilaments and both 
stationary and mobile crosslinkcrs. Both mobile and stationary crosslinkcrs induce filament 
alignment promoting either polar or nematic order. In addition, mobile crosslinkers, such as 
clusters of motor proteins, exchange forces and torques among the filaments and render the 
homogeneous states unstable via filament bundling. We start from a Smoluchowski equation 
for rigid filaments in solutions, where pairwisc crosslink- mediated interactions among the 
filaments yield translational and rotational currents. The large-scale properties of the system 
are described in terms of continuum equations for filament and motor densities, polarization 
and alignment tensor obtained by coarse-graining the Smoluchovski equation. The possible 
homogeneous and inhomogeneous states of the systems are obtained as stable solutions of the 
dynamical equations and are characterized in terms of experimentally accessible parameters. 
We make contact with work by other authors and show that our model allows for an estimate 
of the various parameters in the hydrodynamic equations in terms of physical properties of 
the crosslinkers. 

I. INTRODUCTION 

Soft active systems are a new and exciting class of complex fluids to which energy is continuously 
being supplied by internal or external sources. Biology provides many examples of such systems 
which include cell membranes, biopolymer solutions driven by chemical reactions, collections of 
living cells moving on a substrate, and the cytoskeleton of eukariotic cells [J. The cytoskeleton is a 
complex three-dimensional network of long filamentary proteins (mainly F-actin and microtubules) 
cross- linked by a variety of smaller proteins P, Q] . Among the latter are clusters of active motor 
proteins, such as myosin and kinesin, that transform chemical energy from the hydrolysis of ATP 
(adenosine tri-phosphate) into mechanical work and are capable of "walking" along the filaments, 
mediating the exchange of forces between them 



The self-organization of motor-filament mixtures has been the subject of recent experiments jj, 



|. Specifically, mixtures of microtubules and associated motor clusters have been studied in 
vitro in a confined quasi- two-dimensional geometry Q|. Complex patterns, including asters and 
vortices or spirals have been observed in these in- vitro experiments as a function of motor and ATP 
concentration Q, Q]. The high frequency mechanical response of active filament solutions which 
are dominated by the bending modes of the filaments have also been studied both experimentally 
and theoretically 0, 0, Q|. The study of the properties of these simplified model systems paves 
the way to a better understanding of the formation and stability of more complex structures of 
biological relevance, such as the mitotic spindle formed during cell division Q]. 

There have been a number of recent theoretical studies of the collective dynamics of rigid 
active filaments. First and most microscopic, numerical simulations with detailed modeling of 
the filament-motor coupling have been used to generate patterns similar to those found in experi- 
ments 0,Q|. These approaches have given valuable insights into the problem but are limited to small 
system sizes by computing power. A second very interesting development has been the proposal of 
'mesoscopic' mean-field kinetic equations governing the dynamics of individual filaments where the 



effect of motors was incorporated via a motor-induced relative velocity of pairs of fi 
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. Finally, 



the form of such velocity inferred from general symmetry considerations 
hydrodynamic equations have also been proposed where the large scale dynamics of the mixture 
is described in terms of a few coarse-grained fields whose dynamics is also inferred from symmetry 

considerations, 0, 0, Q, Q, H, Q, Q Q- 

Recently, a connection between the mesoscopic 
and hydrodynamic approaches was established by us by deriving hydro dynamical equations via a 
coarse-graining of the kinetic equations |3j . This was done in the spirit of polymer physics which 
has been successful at predicting macroscopic dynamical behavior of polymer solutions based on 
models of the microscopic dynamics. To make a link with the motor properties we consider a 
simplified model of the motor filament interaction in Appendix EI 

The richness of the phenomena exhibited by the cytoskeleton is illustrated by the ability of its 
constituents to organize in a variety of different structures. In addition different constituents can 
form very similar structures. This leads naturally to the question - how much of the behavior is 
specific and how much is generic? To answer this question it is important to make the connection 
between microscopic models and 'generic' hydrodynamic approaches. 

In this paper we describe a derivation of the hydrodynamic equation for a solution of polar 
filaments and both stationary and mobile crosslinkers. A brief summary of the approach and some 
of the results have been presented earlier [2^]. The filaments are modeled as rigid rods of 
fixed length. Hydrodynamics is obtained by coarse-graining the Smoluchowski equation for rods 
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in solution, coupled via excluded volume and motor-mediated interactions. Small protein clus- 
ters crosslinking the filaments can be grouped in two classes. The first class comprises stationary 
crosslinkers, such as a-actinin, that can induce rotation and alignment of the filaments even in 
the limit of vanishing ATP consumption. Such passive crosslinkers may be polar or nonpolar in 
nature depending on whether they preferentially bind to pairs of filaments of the same polarity or 
their binding rate is independent of the filaments' polarity. They always induce filament alignment 



via a mechanism that has been referred to as "zipping" effect in the literature |2(|. In general we 
expect that most crosslinkers will be polar, although "disordered" motor clusters (i.e, cluster with 
no spatial order in the arrangements of individual motors as in e.g. small myosin clusters) can 
crosslink filaments regardless of their relative polarity. Stationary crosslinkers can lead to the onset 
of the homogeneous nematic and polarized states. The interplay between these two types of order 
is determined by the crosslinkers' polarity. The second class consists of clusters of motor proteins 
crosslinking two filaments, "active crosslinkers". These can also drive the system into nematic and 
polarized states. However, in addition by consuming ATP, the motor heads can "walk" along the 
filaments and mediate the exchange of forces between filaments, inducing filament motion relative 
to the solution (treated here as an inert background). The motor activity depends crucially on the 
ATP consumption rate, which is the driving force that sets up and maintains the nonequilibrium 
state and enters the equation as a chemical potential. Motor activity destabilizes the homogeneous 
states and induces the formation of spatially inhomogeneous structures on mesoscopic scales, remi- 
niscent of those seen in the in vitro experiments. There are two main motor-mediated mechanisms 
for force exchange among the filaments. First, active crosslinkers induce bundling of filaments 
, building up density inhomogeneities. This is the main mechanism responsible for instabilities. 
It is effective only if the rate at which motor clusters step along the filament is inhomogeneous, 
which can be due to crowding and fluctuations in the density of bound motors, or to stalling at 
the polar end. In addition, active crosslinkers sort the filaments according to polarization at a rate 
proportional to the mean motor stepping rate. This mechanism is important in the polarized state, 
where it yields filament advection along the direction of polarization and allows for the onset of 
oscillatory structures. 

The forces and torques exchanged by filaments via the crosslinks are described by considering 
the kinematics of two filaments crosslinked by a single protein cluster that can rotate and translate 
as a rigid object relative to the filaments. The hydrodynamic equations are then obtained by 
suitable coarse-graining of the Smoluchowski equation. This method yields a general form of 
hydrodynamics which incorporates all terms allowed by symmetry, yet it provides a connection 
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between the coarse-grained and the microscopic dynamics. By comparing the equations obtained 
here to those obtained from a microscopic model of the forces exchanged between motors and 
filaments we can relate some of the parameters in the hydrodynamic equations to parameters that 
can be controlled in experiments. 

The hydrodynamic equations are then used to describe the dynamics of the isotropic, nematic 
and polarized solutions. We characterize the possible homogeneous states of the system in term of 
experimentally accessible parameters and discuss the various mechanisms by which motor activity 
can destabilize each homogeneous state. 

In Section [H] we describe the kinetic model of rods crosslinked by small protein clusters and 
set up the formalism of the Smoluchwski equation. The dependence of the crosslinked-induced 
rotational and translational velocities of the filaments on filament orientation and position is ob- 
tained from general symmetry considerations and conservation laws. The details of the kinematics 
of motors and filaments are described in Appendix where a specific microscopic model of the 
coupling is also presented. In Section II I II we obtain the hydrodynamic equations for the system 
by a systematic coarse graining of the Smoluchowski equation. The full form of the hydrodynamic 
equations, including diffusive, excluded volume and active contributions, is given in Appendix O 
The nonlinear hydrodynamic equations are solved in Section IIVI to obtain the possible homoge- 
neous steady states of the system. "Phase diagrams" are constructed in terms of the filament 
and crosslinkers densities identifying the isotropic, nematic and polarized states. The nonlinear 
hydrodynamic equations for each homogeneous state are presented in Section where the sta- 
bility of each state is also studied. All homogeneous states become unstable at high filament and 
crosslinkers densities via filament bundling. The interplay of bundling and diffusion promotes the 
onset of stable spatial structures on mesoscopic scales. Finally, we conclude with a discussion of 
open questions and a comparison with related work. 

II. THE MODEL: SMOLUCHOWSKI EQUATION FOR MOTOR-FILAMENTS 

SOLUTIONS 

We model the system as a collection of thin rods of fixed length I and diameter b « I crosslinked 
by small protein clusters (of linear size ~ b) that can exchange torques and forces between the 
filaments. Filaments and crosslinkers move through a solvent which is assumed inert. The solution 
forms a quasi-two-dimensional film, of thickness much smaller than the length of the filaments. The 
dynamics of both filaments and crosslinkers is overdamped. This is a good model for a quiescent 
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solution with no externally imposed flow nor net flow generated by motor activity. We are interested 
in describing the filament dynamics on time scales large compared to the characteristic times for 
binding and unbinding of the crosslinkers so that we can treat a constant fraction of them as bound. 
The dynamics of crosslinkers binding and unbinding was considered for instance in Ref. [2| and 
it was found that varying the rates of binding and unbinding of motor clusters did not affect the 
nature of the nonequilibrium steady states of the active solution. The temperature of the system 
is taken to be constant and the effect of thermal fluctuations is not considered explicitly. We 
assume, however, that the stochastic nature of the crosslinkers dynamics, as well as other sources 
of noise in the systems, can be incorporated in an effective temperature T a that may differ from 
the actual temperature of the solution Finally, although the kinetic model described below 

applies to a solution with a low concentration of filaments, the structure of the continuum equations 
obtained upon coarse-graining the kinetic model is general and not restricted to low density. On the 
other hand, the quantitative estimates obtained for the various parameters in the hydrodynamic 
equations are for a low density of filaments and crosslinkers. 

The dynamics of the concentration c(r, u, t) of filaments with center of mass at r and orientation 
u at time t is governed by the Smoluchowski equation [271 . |28| |. which describes conservation of 
the number of filaments, 

d t c= -V- J c -K- J c , (2.1) 

where TZ = u x d-a is the rotation operator. The translational current density, J c , and rotational 
current density, J7" c , are given by 

J ci = -D l3 V jC - — VjV cx + 4 , (2.2) 
J cl = -D r K iC - -^-cTZiV cx + j£ , (2.3) 

KBT a 

where Dij = DuUiiij + D±(5ij — UiUj) is the translational diffusion tensor and D r is the rotational 
diffusion rate. For a low-density solution of long, thin rods D± = Dn/2 = D/2, where D = 
kBT a ln(l/b)/(27TT]l), with r/ the solvent viscosity, and D r = 6D/1 2 . The potential V cx incorporates 
excluded volume effects which give rise to the nematic transition in a solution of hard rods. It can 
be written by generalizing the Onsager interaction to inhomogeneous systems as ksT a times the 
probability of finding another rod within the interaction area of a given rod (see Figure ^) . In two 
dimensions this gives 

y cx (ri,iti) = k B T a / dr 2 \ du 2 c(r 2 ,u 2 ,t) \u± x u 2 \ / 5(r 1 + tiisi - r 2 - u 2 s 2 ) 

J J J S±S2 




FIG. 1: The geometry of overlap between two interacting filaments of length I cross-linked by an active 
cluster. The cross-link is a distance si, (S2) from the centre of mass of filament 1(2). The distance between 
centres £ = T2 — r% = sxiii — S2«2- 



k B T a / du 2 / |«i x u 2 \ c(ri + £,u 2 ,t) , 
J J S1S2 



(2.4) 



where Sj, with —1/2 < s» < Z/2, parametrizes the position along the length of the i-th filament, for 
z = 1,2, and J s ... = /_j/ 2 dsi.... The 5-function ensures that the filaments be within each other's 
interaction volume, i.e., in the thin rod limit b « I considered here, have a point of contact. 
The factor |tti x 62 1 represents the excluded area of two thin filaments of orientation ui and U2 
touching at one point. In the second equality we let £ = v 2 — ri = U\S\ — U2S2 j^ . 

The translational and rotational active current of filaments with center of mass at ri and 
orientation along ui are written as 



Jc(ri,tti) = b 2 I I \ui x u 2 \m(ri + uisi)vi(si,s 2 ,ui,u 2 )c(ri,ui,t)c(ri + £,u 2 ,t)(2.5) 
J^{r x ,ui) = b 2 / / \ui x u 2 \m(ri + uis^u^si, s 2 ,iii,u 2 )c(r 1 ,u 1 ,t)c(r 1 +£,u 2 ,tj2.6) 

JU2 J S\S2 

where mix) is the density of bound crosslinkers, evaluated at the point of attachment to the 
filaments. Also, vi = ri is the translational velocity that the center of mass of filament 1 acquires 
due to the crosslinker-mediated interaction with filament 2, when the centers of mass of the two 
filaments are separated by £. Similarly, u}% X u\ = u\ is the crosslinker-induced velocity of rotation 
about the center of mass (see Figure 
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Small clusters of motor proteins crosslinking two filaments can be stationary or mobile (active). 
In general the density of bound crosslinkers can be written as m = m s + m a , with m s and m a the 
densities of stationary and mobile crosslinks, respectively. Mobile crosslinks are clusters of motor 
proteins that can diffuse and be convected along the filaments. The mobile crosslinker density 
obeys a diffusion-convection equation given by 

d t m a = D m V 2 m a - V • J m , (2.7) 

where 

and u(s) is the speed at which a motor cluster steps on a filament at position s. The mean 
value «o = f s u ( s )/l °f the stepping rate is uq ~ ciRatp, where a is the step size and Ratp 
is the ATP consumption rate. For typical motor clusters (kinesins on microtubules or myosins 



/ / u u(s) c(r, u, t)m a (r + us, t) , (2i 

J u J s 



on F-actin) no ~ nm/msec P). As shown in Ref. |29|, spatial inhomogeneities in the motor 
stepping rate u(s), especially the stalling of motors at the polar end, are crucial for driving filament 
bundling and pattern formation. Such inhomogeneities may arise from motor crowding or from 
large fluctuations in the concentration of ATP under condition of near depletion. Very recent 
experiments in purified actin-myosin II solutions have indeed suggested that the motor-driven 
formation of spatially inhomogeneous patterns, such as asters and vortices, may be associated with 
strong inhomogeneities in motor activity j^ . 

The translational and rotational velocities of the filaments induced by crosslinkers are written 
in a general form that is consistent with translational and rotational invariance. We consider a pair 
of filaments cross-linked by a single protein cluster. As seen below, all crosslinkers can exchange 
torque among the filament and induce filament alignment or "zipping". Mobile crosslinkers that 
consume ATP to step along filaments can also exchange forces and induce translational motion of 
the filaments. In general the rotational and translational dynamics induced by the crosslinkers is 
coupled. 

It is convenient to introduce the relative velocity and net velocity of the filament pair as 

V = Vi - v 2 , 

V = , (2.9) 

with vi ; 2 = V ± v/2. A general form of the relative linear velocity v and angular velocity u = 
u\ — u>2 of the filament pair consistent with symmetries and conservation laws is 

^m (+ m^. ai) , (2 . 10 , 
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u = 47(61) uixu 2 , (2.11) 

where £ = «isi — U2S2 is the separation of the filaments' centers of mass, and a, (3 and 7 depend 
on the relative orientation of the two filaments through the angle 6 = cos~ 1 (iii • U2). The angular 
dependence of a, (3 and 7 arises both from the kinematics of the crosslinker-mediated filament 
interaction, as well as from the dependence of the probability that a protein cluster binds two 
filaments on the angle between the filaments at contact. 

It is instructive to rewrite the relative velocity v in terms of two orthogonal vectors as 

5(0), WA A , \P{0) ale). .-1,* A , , 

v = -^(si - s 2 )(u 1 + u 2 ) + [Zy - -^-(8! + s 2 )\ (u 2 - u x ) . (2.12) 

The physical meaning of Eq. (|2.12j) can be understood by considering a specific microscopic model of 
the motor-filament coupling, such as the one described in Appendix^ In this model the kinematics 
of two filaments coupled by a motor cluster is described explicitly in terms of the rate u(s) at which 
the cluster steps along the filament and the torsional stiffness k of the cluster A comparison 
of Eq. (|2.12|) with Eq. (|Allj) assuming a linear dependence of u(s) on s as u(s) ~ no — u's, with 
v! = —du/ds, shows that in the microscopic model a and (3 are independent of the angle 0, with 
(3 = uq and a = 2lv! . In general we can identify (3 with the mean rate at which a motor cluster 
steps along a filament, i.e., (3 ~ (l/l) f u(s), while a is controlled by spatial variation in the 
stepping rate, with a ~ 21 max \du/ds\. It is then apparent that the first term on the right hand 
side of Eq. (|2.1Ujl arises from variation in motor activity along the filament, such as the stalling 
of motors before detaching upon reaching a particular point on the filament. It is proportional to 
the separation £ of the filaments' centers of mass and vanishes when these coincide. The angular 
dependence of a is chosen so that this contribution to the relative velocity is largest when filaments 



are parallel. The second term in Eq. 1)2. lOj) . proportional to [3, vanishes for aligned filaments and 
drives the separation or sorting of anti-aligned pairs. 

For small angles, we can write the functions a and (3 and 7 in the form of expansions in powers 
of Ui ■ U2 as 

ol(0) ~ a + &i(ui ■ u 2 ) , (2.13) 
(3(9) - Po + f3i(ui ■ u 2 ) , (2.14) 

where all coefficients are defined positive. It can be shown that within the approximation used 
below, where we only consider the first three moments of the filament concentration, no new terms 
are obtained in the continuum equations for such moments when terms of higher order in Ui U2 are 
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included in Eqs. Q2.13|) and (|2.14[) . Contributions of higher order in the angle between the filaments 
only affect the numerical coefficients of the various parameters in the continuum equations. 

The rotational parameter 7 can be estimated by describing the crosslinker as a torsional spring 
of constant k, as shown in Appendix where we find that the rotational rate induced by a single 
crosslinker does not depend on 9. We estimate 7 ~ D r n/kBT a . In this model we assume that the 
motor cluster always binds on the side of the smaller angle between the filaments, as shown in Fig. 
12 We distinguish between polar clusters that bind preferentially to filaments of the same polarity 
(Fig. Efa)) and nonpolar clusters that bind to filaments regardless of their relative polarity (Fig. 
|2fb)). The probability for such two classes of protein clusters to bind to filaments will in general 
depend on the angle 9 between the filaments, yielding an angular dependence of the effective rate 
7(0). Again, to lowest order in u\ ■ U2 we write 

^ IP + 1np{u\ ■ u 2 ) , (2.15) 

The term proportional to 7p favors rotations that align filaments of the same polarity and describes 
polar clusters 0,B|, which are in general expected to be active crosslinks in the presence of ATP. 
The term proportional to 'Jnp favors rotation in the direction of angles 9 < ir, regardless of the 
relative polarity of the two filaments. It describes non-polar clusters which bind to filament pairs of 
any orientation Passive cross- linkers (such as a-actinin on F-actin which play a crucial role in 
the rheology of actin gels) j^l can be either polar or nonpolar. Polar clusters (jp 7^ 0) where not 
considered in earlier work by two of us j^ . but are crucial for the formation of a polarized phase 
(see also Ref. j3|). Both 7p and 'Jnp will increase with increasing binding rate of the clusters to 
the filament. It is interesting to speculate that the kinesin constructs in the experiments by Nedelec 
et al. 0, Q] are polar clusters, while the disordered myosin II clusters studied by Humphreys et 
al. P| may be apolar in nature. We can also imagine that if the binding/unbinding of the motor 
clusters does not require ATP, these terms, unlike the active contributions to the translational 
currents, would be independent of the ATP hydrolysis rate. 

To determine the net linear velocity V and rotational velocity fi = (u\ + ^2)j1 we note that 
the third law and momentum conservation require that the net force and torque due to a motor 
cluster on a pair of filaments vanishes in the absence of external forces. This yields 

Cij(ui)v!j + (ij(iL2)v2j = , (2.16) 
CrOJl + Cr<^2 = , (2.17) 

where Ciji^ 1 ) = C\\UiUj + C±($ij ~ ^iUj), with 0| = ksT a /D^ and (j_ = ksT a /D^, is the friction 
tensor of a long thin rod and Q T = k£,T a /D r is the rotational friction. The vanishing of the net 
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(a) 



(b) 



FIG. 2: Polar and non-polar clusters interacting with polar filaments. Assuming that clusters always bind to 
the smallest angle, polar clusters ('Jp/^np 3> 1) bind only to filaments in configuration (a) while non-polar 
clusters (jp/jnp <C 1) bind to both configurations equally. 

torque on the pair clearly requires u 2 = — uf±, i.e., there is no net rotational velocity. The net 
velocity V of the pair is generally nonzero and is given by the solution of Eq. (|2.1fi|) . or 

1, 



[Cij(«l) +Cij(u 2 )]Vj 
Its explicit form is given by 

V = A(ui + u 2 ) + B(u 2 - iii 

with 

.4 
B 



a 1 — u\ ■ u 2 ' 



4 1 — au\ • u 2 
a 1 + u\ ■ u 2 a{9) 



:o)-^(s 1 + S2 ) 



21 

Oi - s 2 ) 



(2.18) 

(2.19) 

(2.20) 
(2.21) 



4 1 + ou\ ■ u 2 21 

where a = ((± — C||)/(Ci + Cll) > 0- F° r l° n § thin rods (j_ = 2^|| and a = 1/3. 

Equations ()2.12j) and (|2.19|) display explicitly the even and odd symmetry of V and v, respec- 
tively, under filament exchange. Note that the net velocity V vanishes for isotropic bodies, i.e., 
when <ji = C± { a = 0)- 



III. CONTINUUM EQUATIONS 

Our goal here is to obtain a set of coarse-grained equations to describe the macroscopic dynamics 
of active filament mixtures on scales large compared to the filaments' length, I and on timescales 
long compared to the typical binding times of the cross-linkers. 



This level of description is valid when the macroscopic quantities describing the so 
hibit spatial variations on length scale much greater than the length of the filaments 



ution ex- 



2J. The 
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macroscopic quantities we choose to study are the local filament density, p(r,t), the local filament 
polarization, p(r,i), and the alignment tensor, Sij(r,t), a second rank symmetric tensor which 
measures the local orientational order in a nematic state. These fields are associated with either 
conservation laws (the density) or possible broken continuous symmetries (p, SV,) and therefore 
control the hydrodynamic modes of the system. They can be defined as the first three moments of 
the distribution c(r, u,t) j^j], 

P(r,t) = 
T(r,t) = 
Qij(r,t) 

with Qij(u) = UiUj — Hydrodynamic equations for these coarse-grained densities can be 

obtained by writing an exact moment expansion for c(r, u, t) (see Appendix[B|) and truncating this 
expansion at the third moment. To derive the continuum equations we assume that all quantities 
are slowly varying on the scales of interest and expand the concentration of filaments c(ri +£,1*2) 
and the crosslinker density m(rx + U1S1) in the expressions for the active currents near their values 
at ri as 

c(ri + £,u 2 ) = c(n,u 2 ) +tidi i c(r 1 ,u 2 ) + -&&<9iidijc(ri,«2) + 0(V 3 ) , (3.2) 

m(r 1 + u\si) = m(ri) + uiiSi3um(ri) + ^uiiUijs\diidijm(r 1 ) + C(V 3 ) . (3.3) 

When the expansions (|3.2|) and 1)3.3 J) are inserted in Eqs. 1)2.5)1 and 1)2. 6)1 . the integration over s± and 
S2 can be carried out term by term. An analogous expansion is used to approximately evaluate 
the excluded volume interaction, as well as in the equation for the motor concentration. Some 
details of the derivation of the hydrodynamic equations for the motor density, filament density, 
polarization and alignment tensor are given in Appendix \d 

For simplicity, here we consider the nonlinear continuum equations retaining only terms up to 
second order in the gradients. While the analysis of the linear stability of homogeneous states with 
terms up to fourth order in the gradients does introduce a new length scale (see Appendix [DJ), 
the simplified equations are instructive and capable of describing much of the physics. The motor 
density obeys a simple diffusion equation given by 

d t m a = D m V 2 m a - V • (m T) , (3.4) 



/ c( r,u,t) , 

J u 

e p(r,t) p(r,t) = u c(r,u,t) , 

J u 

= p(r,t) Sij(r,t) = / Qij(u) c(r,u,t) , 



(3.1) 
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where the second term describes convection of the motors along the filaments ^?|]. The equations 
for the filament density, polarization and alignment tensor are 

dtp = -diJi , (3.5) 

dt(pPi) = -djJij - Ri , (3.6) 

dt(pSij) = -d k J ijk - Rij , (3.7) 

where the currents are given by 

Ji{r,t) = I Jd(u,r,t) , (3.8) 

J u 

Jij( r ,t) = UiJ cj (u,r,t) , (3.9) 

J u 

Jijk(r,t) = / QijJ ck (u,r,t) . (3.10) 

J it 

The rotational current does not contribute to the density equation, but it yields the source terms 
Ri and Rij in the equations for the polarization and alignment tensor. These are given by 

Ri(r, t)= I UiK- J c (u, r, t) , (3.11) 
Rij{r,t) = I Qij n-J c (u,r,t) . (3.12) 

J u 

The explicit form of the translational (|3.8I3.1U|) and rotational 1)3.1113.12(1 currents is given in Ap- 
pendix 10 The equation for the density p has the form of a continuity equation, as required by 
filament number conservation. The local polarization p and the alignment tensor S$j define the 
order parameters needed to characterize the ordered states of the system and are not conserved 
variables. Each ordered state discussed below will, however, be characterized by a broken orienta- 
tional symmetry and a corresponding broken symmetry variable (a unit vector along the direction 
of broken symmetry) whose fluctuations are infinitely long lived at large wavelength, as required 
for hydrodynamic modes. 
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IV. HOMOGENEOUS STATES 

We begin by identifying the possible homogeneous steady states of the system. In this case 
all contributions to the dynamical equations for the filament solution come from the rotational 
currents. Both the motor and filament densities have constant values m = mb 2 and p = p$. The 
equation for the polarization and the alignment tensor are given by 

d t pi = -[D r - j P mp ]pi + 4D r p /p N + (j NP - 2<y P )rhp SijPj , (4.1) 



dtSij 



4D r (l - Po/pn) ~ iNprhpo Sy + 2j P mp (piPj - ^p 2 ) , (4.2) 



where all filament densities are measured in units of I 2 , and pn = Sir/2. The motor-induced 
rotational rates -fp and 7_/vp have dimensions of frequency and represent the effect of polar and 
nonpolar motor clusters, respectively. For simplicity we denote by rh the total dimensionless 
density of crosslinkers, without distinguishing between stationary and active protein clusters. One 
can imagine situations, however, where jp will in general be proportional to the ATP consumption 
rate, but the nonpolar coupling ^np will be only weakly affected by ATP concentration. In the 
following all lengths are measured in units of the filament length I and times are measured in units 
of D' 1 . 

There are three possible homogeneous stationary states for the system, obtained by solving 
Eqs. (|4~Tj) and (|4~2j) with d t pi = and d t Sij = 0. These are: 

isotropic state (I) : pi = Sij = , 
nematic state (N) : Pi = Sij ^ , 
polarized state (P) : Pi ^ Sij ^ . 

At low density the only solution is pi = and Sij = and the system is isotropic (I). The homo- 
geneous isotropic state can become unstable at high filament and/or motor density, as described 
below. 

To discuss the instabilities it is convenient to rewrite Eqs. (|4.1(l and (|4.2j) in a more compact 
form as 

d t pi = -am + b 1 poS ij pj , (4.3) 
dtSij = -a 2 Sij + b 2 po(piPj - ^p 2 ^ ■ (4.4) 
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The coefficients a\, b\, 02, and 62 are given by 

ai = 1 - ihjppo/Dr , (4.5) 

a 2 = 4[1 - Po/pn ~ 7Nprhp /(4D r )] , (4.6) 

61 = 4[p^ + ( 77V p - 2 lP )rh/{AD r )} , (4.7) 

6 2 = 2-yprh/Dr . (4.8) 

In the absence of crosslinkers (7P = ^np = 0), no homogeneous polarized state with a nonzero 
mean value of p is obtained. There is, however, a transition at the density pn = 3n/2 from an 
isotropic state with Sij = to a nematic state with Sij = SofaiUj — \8ij) , with n a unit vector along 
the direction of broken symmetry. The transition here is identified with the change in sign of the 
coefficient a 2 of Sij on the right hand side of Eq. 1)4. 4j) . A negative value of 02 that controls the decay 
rate of Sij signals an instability of the isotropic homogeneous state. This occurs when excluded 
volume effects dominate at po = pjy- The homogeneous state is isotropic for po < pn and nematic 
for po > Pn- A mean- field description of such a transition, which is continuous in 2d (but first order 
in 3d), requires that one incorporates cubic terms in Sij in the equation for the alignment tensor. 



Adding a term -c 2 plS k iS k iSij to Eq. (jOJ) we obtain So = ^ \J ' ~2a 2 / c 2 = ±^-8(1 - Pq/pn)/c2- 
If 7p = 0, but 7jvp 7^ 0, there is again no stable polarized state. The presence of a concentration 
of nonpolar crosslinkers does, however, renormalize the isotropic-nematic (IN) transition, which 
occurs at a lower filament density given by 

p IN (m) = — — — ; . (4.9) 

1 + m-f N pp N /(AD r ) 

The presence on nonpolar crosslinks favors filament alignment and shifts pin downward, as shown in 
Fig. |21 A qualitatively similar result has been obtained in numerical simulation of a two-dimensional 



system of rigid filaments interacting with motor proteins grafted to a substrate |2J] . In this case the 
motors promote alignment by exerting longitudinal forces on the filaments. The amount of nematic 
order So is also enhanced by motor activity, with Sq = — yf— 202/02 = — a/— 8(1 — po/ pin{™))/c2- 
If 7p is finite, the system can order in both polarized and nematic homogeneous states. The 
homogeneous isotropic state can become unstable in two ways. As in the case 7p = 0, a change 
in sign of the coefficient 02 signals the transition to a nematic (N) state at the density pi^{m) 
given in Eq. (|4.9|) , In addition, the isotropic state can become linearly unstable via the growth of 
polarization fluctuations in any arbitrary direction. This occurs above a second critical filament 
density, 

pip{m) = , (4.10) 

7 pm 
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defined by the change in sign of the coefficient a\ controlling the decay of polarization fluctuations 
in Eq. ([4.3)1 . For po > Pip(fh) the homogeneous state is polarized (P), with p / 0. The alignment 
tensor also has a nonzero mean value in the polarized state as it is slaved to the polarization. 
The location of the boundaries between the various homogeneous states is controlled by the ratio 
g = 7p/7tvp that measures the polarity of motor clusters. One can identify two scenarios depending 
on the value of g. 

I) For g < 1/4, the density pip is always larger than pm and a region of nematic phase exists 
for all values of to. At sufficiently high filament and motor densities, the nematic state becomes 
unstable. To see this, we linearize Eqs. (|4.3j) and (|4.4|) by letting = Sfj + 5Sij and 5pi = pi. 
Fluctuations in the alignment tensor are uniformly stable for a 2 < 0, but polarization fluctuations 
along the direction of broken symmetry become unstable for a\ < pob\So/2, i.e., above a critical 
density 



D r 

Pnp = — 

mgjp 



l + -^(l-,ll + 2C2R{1 - R) 



(4.11) 



C2 R v y b( 

where R = pin/ Pip- The polarized state at po > p^p has 

P°i=PoPi, (4.12) 

Sfj = SpiPiPj - 5^/2) , (4.13) 



with p a unit vector in the direction of broken symmetry and 

2d! 



2 2aia 2 
Po 



Pohb 2 L \hpoSo 



(4.14) 



Sp = sJl-f^pl = 2^- . (4.15) 
V 2a x a 2 poh 

The "phase diagram" for g < 1/4 is shown in Fig. 

II) When g > 1/4, the boundaries for the I-N and the N-P transitions cross at 

_ pnDt/jp 

mx ~ 1 - l/(4<7) ' { ' 

where pin = pip = p^p and the phase diagram has the topology shown in Fig. ^ For to > rh x the 
polarity of motor clusters renders the nematic state unstable at all densities larger than /Ojat(to) 
and the system goes directly from the I to the P state at pip, without an intervening N state. 
At the onset of the polarized state the alignment tensor is again slaved to the polarization field, 
Sij = TTPo (PiPj ~ \^ijP 2 ) > an d P = Pop- The value of po is determined by cubic terms in Eq. 
()4.3|) not included here. 
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NP 



~ m 



FIG. 3: The homogeneous phase diagram for g < 1/4. For all values of m a region of nematic phase exists 
between the isotropic and polarized phases ("fp/D r = 1, g = 1/10 and c 2 = 50). 



P 




FIG. 4: The phase diagram for g > 1/4. For m > rh x , where pjjv and p/p intersect, no N state exists and 
the system goes directly from the I to the P state (jp/D r = 1, g = 1 and c 2 = 50). 



Finally, we note that if 7jvp = 0, with 7p 7^ (i.e., g — ► 00), the I-N transition is independent 
of motor density and always occurs at po = Pat- The motor density where pin = pip reduces to 
rh x = p N D r /j P . 

Estimates of the various parameters can be obtained using a microscopic model of the motor- 
filament interaction of the type described in Appendix Using parameter values appropriate for 
kinesin (k ~ 10~ 22 nm/rad |2() we estimate 7p ~ 7^p ~ i^/Cr = KD r /(kp.T a ) ~ 10~ 1 sec~ 1 , where 



17 



we used the value D r ~ 10 -2 sec -1 appropriate for long thin rods in an aqueous solution [351 ] and 
T a ~ 300 K. Using / ~ 10pm, b ~ 10 nm, the value fh x above which no nematic state exist is 
found to correspond to a three-dimensional crosslinker density of about 0.5 — lpM for g = 1 and a 
sample thickness of order lpm. This value is of order of the motor densities used in experiments 
on purified microtubule- kinesin mictures such as those of Ref. Q|, suggesting that the filament 
solution in this experiments is always in the polarized state. 

V. DYNAMICS OF INHOMOGENEOUS STATES 

Experiments on motor /filament mixtures have shown that uniform states are often unstable to 
the formation of complex spatial structures. These instabilities arise from the growth of spatial 
fluctuations in the hydrodynamic fields In particular, the rate of motor-induced 

filament bundlingcan exceed that of filament diffusion yielding the unstable growth of density 



inhomogeneities 
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MUM- 

States with spatially varying orientational order, where 




33]. 



the filaments spontaneously arrange in vortex and aster structures, are also possible |17l . |23j, 
To understand the different nature of the instability from each homogeneous state, we now examine 
the dynamics of spatially- varying fluctuations in the hydrodynamic fields in each of the stationary 
homogeneous states of the system. The hydrodynamic fields are those with characteristic decay 
times that exceed any microscopic relaxation time and become infinitely long at long wavelengths. 
We find that the low frequency hydrodynamic modes of this active system are determined by 
fluctuations in the conserved densities and in variables associated with broken symmetries. A 
change in sign in the decay rate of these modes signals an instability of the macroscopic state of 
interest. For simplicity we only discuss here the case of constant motor density. In addition in this 
section we let ao = ao/48, a\ = di/48 and assume ao = a± = a, Po = Pi = P- This approximation 
is justified by the estimate for these motor-induced parameters obtained in Appendix [3 using a 
microscopic model for the the motor-mediated filament interaction [^J. We consider separately 
the spatially varying hydrodynamic modes in each of the homogeneous states: isotropic, nematic 
and polarized phases. 

A. Isotropic state 



The isotropic homogeneous state has p = po, po = and Sfj = 0. The only hydrodynamic 
variable in this state is the density of filaments. The polarization vector p and nematic orientation 
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tensor Sij are not hydrodynamic modes and therefore relax to zero on microscopic timescales. 
However, due to dynamical constraints such as entanglements the relaxation of p and Sa can 
become sufficiently slow to yield finite lifetime and finite wavelength inhomogeneities 2J, |3j| . 
The nonlinear equation for the density is given by 

3D 

d t p = —V-(l + v p)Vp-aV-{pVp), (5.1) 

where vq = 2/ir. The active current proportional to a has an effect opposite to that of thermal 
diffusion as it tends to build up density inhomogeneities. As we will see below this term drives 
filament bundling and is the main pattern-forming mechanism in each of the homogeneous states. 

1. Linear Stability 

To examine the stability of the isotropic state we consider the dynamics of fluctuations 5p(r) = 
p(r) — po of the density about its steady-state value, po to linear order in 5p. By expanding the 
fluctuations in Fourier space, 

Mr) = S kPk e ik - r , (5.2) 

the linearized equation for the Fourier amplitudes is given by 

r3 i 
dtPk = -k 2 -D(l + p v ) - ampo pk • (5.3) 

The relaxation of density fluctuations is governed by a diffusive mode of frequency 

r3 i 
z p (k) = -k 2 -D(l + p t>o) - amp . (5.4) 

The isotropic state becomes unstable against the growth of density fluctuations if z p (k) > 0, or 
a > a c , with 

3D(1 + po^o) 

° c = 4mp • (5 ' 5) 
Conversely, the homogeneous isotropic state becomes unstable for filament densities larger than a 
critical value 

j 3D 3D . . 

p l B = . (5.6) 

Arha — 3Dvq Ama 



I, 29]. 

description of the bundling instability requires that one incorporates terms up to 4 th order in the 
gradient expansion of the hydrodynamic fields j^fj]. The 4 th order terms introduce a new length 
scale above which the homogeneous state becomes again stable, as shown in Appendix [Dl The 
onset of the instability is, however, controlled entirely by the quadratic terms considered here. 
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B. Nematic state 



The homogeneous nematic state is characterized by p = po, po = 0, and Sij = So(nirij 



where n is a unit vector in the direction of broken symmetry, known as the director field. For 
concreteness we choose n = y. The hydrodynamic fields of such an overdamped nematic liquid 
crystal are the density and the director. The symmetry of the nematic state requires that the 
hydrodynamic equations incorporate the symmetry n — > — n. The magnitude So of the alignment 
tensor is not a hydrodynamic field and will be assumed constant below. For simplicity we also 
neglect excluded volume corrections. The nonlinear hydrodynamic equations for filament density 
and director field are then given by 



d t p = ^(^1 - y^V 2 p - ^(1 - S'o)amV 2 / o 2 + Y<9i (D — 4amp)r 



+ Y 9i 



(D - -aihp)pdj(ninj) 



(5.7) 



dm 



7D 1 / S 

2^ + 8l 1 -y )wpm 



V 2 (pn J ) + [— + hi - So^Nprh] djV ■ (pn) 
L p 8 J 



+ 
+ 

+ 



— + hi + 2S )'j N prh n k nidid k (pnj 
p 8 J 

1" o 1 + tt )lNpm 



^(l-^-\ + hl-^--^ 



p \ 

D 1 

7 + 8 



4V 
1 



1 



So 



2S 



iNpm 



n k d k djp 



So 

T 



iNpm 



pd k (n k m)dinj 

^^am|^— n fe 9 fc 5jp 2 + 5S pd k (n k ni)dinj + 2Sopni(din j )d k n k 

7 9 I") 

(5 - 3So)d k n j + -djn k + -S jk dini + 2(2 + 3S )n k nidinj d k pj 



pn k 



(d k ni)djni^ 



(5.8) 



where 8f- = 5ij — riirij projects in the direction transverse to n. 

In the density equation, Eq. (|5.7|) . activity plays the same role as in the isotropic state, with 
bundling (~ a) opposing diffusion and eventually driving the instability of the homogeneous state, 
as described below. The first and second terms on the right hand side of Eq. (|5.8|) for the nematic 
director are the elastic restoring forces associated with bend and splay deformations, respectively. 
These elastic constants are softened by filament bundling, while motor-induced alignment (~ jnp) 
tends to stabilize the homogeneous nematic state. A solution with p = constant requires vanishing 
of splay, i.e., V • n = 0. In this case the director equation reduces to 



dm 



7D 1 / S 
— + -[l--)lNPmp 



So 

+ ^lNpmpn k (d k ni)[djni + dirij] 



1 



V z rij+ D + -(l + 2Soh N pmp 



n k nidid k nj 



10S 



ampn k (d k ni)(dirij)} 



(5.9) 
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In this case bundling does not play any role to linear order. 

To discuss the stability of the homogeneous nematic state, we consider the dynamics of spatially 
varying fluctuations of the hydrodynamic fields about their mean values, by letting 



p(r) = p + 5p(r) , 
n(r) = y + 5n_|_(r) 



(5.10) 
(5.11) 



To lowest order in the fluctuations 5n± is perpendicular to y, i.e., in the two-dimensional geometry 
considered here, 5n± = 5n x x. The linearized equation for the Fourier amplitude of density and 
director fluctuations for Sq = 1 are given by 



dtnk 



dtPk 



1 



1 r 

2 



Dk 2 + (D- 4amp )kl 



D 



1 



1 



arhpo) pokxkyTik 



(5.12) 



3D + —j NP rhp \k + -jNpmp k y n k - -[D - -amp )k x k y — . (5.13) 



Pk 



24 ri 7 4 M,J ru y l ~ 8\ 3 ■ -J " po 

The hydrodynamic modes in the nematic state describe the coupled decay of density and director 
fluctuations. They are always diffusive and are given by 



z±(k,(j>) 



(5.14) 



where 



V±(4>) = ^{^D + —jNpfhpo + (D + -jNprhpo - 4amp ) cos 2 </>} 

=F^{ \ D + j^lNprhpo - (D - ^jNpmpo - 4amp ) cos 2 < 

/ 4 w 8 \ o 2 1 1/2 
+2yD — —ampoj yD — —ampoj sin ^>cos (f> j , 



(5.15) 



and 4> is the angle between the wavevector k and the direction of the broken symmetry (y). To 
gain some insight in the angular dependence of the modes it is useful to consider the behavior for 
special directions of the wavector. For wavevectors k along the y direction (eft = 0), density and 
director fluctuations decouple and we obtain 



Zp(ky) 



D — 2ampo 



k 2 
K y ' 



1 

"4 L 



7 



3D + — jNprhpo 



k: 



(5.16) 
(5.17) 



In this case the director fluctuations are always stable, while the density fluctuations become 
unstable for filament densities above a critical value p^, given by 



D 

2rao 



(5.18) 
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FIG. 5: (color online) The critical filament density Pg((j>) where the homogeneous nematic state becomes 
linearly unstable is shown as a function of the angle (ft between the wavevector k and the direction of broken 
symmetry for D = 1, jnp = 1, fh = 1 and two values of a. At (ft — the critical density is given by 
Eq. I|5.18(l . When either a or m are increased, the density Pb(4>) shifts to lower values at all angles and the 
region of stability of the homogeneous nematic state shrinks. 

For k along x (0 = tt/2) density and director fluctuations again decouple, but both eigenvalues are 
always negative, with 



The homogeneous nematic state is linearly stable for all parameter values against long-wavelength 
fluctuations that only exhibit spatial variation in the direction normal to that of the mean filament 
orientation. 

In general the critical filament density p^(cft) above which the homogeneous nematic state is 
unstable has a complicated angular dependence. It increases with (ft and it diverges for (ft — > 7r/2, 
where the homogeneous state is linearly stable for all filament density. Its angular dependence is 
shown in Fig. El for a few values of parameters. 

For all angles the instability is controlled by the bundling rate, a, while the rotational rate jnp 
always tends to stabilize the homogeneous nematic state. 




(5.19) 



(5.20) 
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C. Polarized homogeneous state 

The homogeneous polarized state is characterized by p = po, po = Pop, and Sfj = Sp(piPj — 
2<%)> with p a unit vector pointing along the direction of broken symmetry. The hydrodynamic 
fields are the filament density and the unit vector p. The magnitude po of the polarization is not 
a hydrodynamic field and will be assumed constant in the following. In the polarized state the 
alignment tensor is slaved to the polarization and does not yield any additional hydrodynamic 
field. Assuming for simplicity po = 1 and neglecting excluded volume corrections, the nonlinear 
hydrodynamic equations for filament density and polarization direction are given by 



and 



where 



dtp ~ ^/3mV • (p 2 P ) 
3d 



3 3 1 

-DV 2 p — -amV 2 p 2 — -amdidj(p 2 piPj) 



txi 
3d 



Pi 



13m 
l!6~ 

am 
~lp~ 

am T 



5D 



+ 



1 

96 

7pm 



D 



24 J ' xrrj ' ' 4p 
$Tj d k{p dj(ppk) + dk(ppj) + 5 jk V ■ (pp) | 

2d k {pp k djp) + d k {ppjd k p) + dj(pp ■ Vp) 



(5.21) 



(5.22) 



<^(p) = $ij -PiPj 



(5.23) 



projects in the direction transverse to p. The usual elastic constants K\ and K3 for splay and 
bend deformations, respectively, can be identified as K% ~ 5D/8 and K\ — K3 ~ D/A. 

The first term on the right hand side of the density equation (Eq. (|5.21|) ) is simply filament 
diffusion. The second term proportional to a opposes diffusion and describes the effect of filament 
bundling. Finally, the last term describes higher order nonequilibrium couplings between density 
and polarization. 

The broken directional symmetry of the polarized state yields an effective drift velocity ~ j3mpp 
describing filament advection along the direction of polarization. This leads to convective-type 
terms on the right hand side of both the density and polarization equations. These are true 
nonequilibrum terms that cannot be obtained from derivatives of a free energy. They arise because, 
due to the anisotropy of rod diffusion, a motor cluster cross-linking two filaments can yield a net 
velocity of the pair, even in the absence of net forces, as shown in AppendixEI This term is absent 
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in descriptions of the hydrodynamics of active polymer solutions and gels close to equilibrium 
proposed on phenomenological grounds on the basis of symmetry argument 
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It is 

therefore a far from equilibrium contribution to active filament dynamics. Kruse et al. |21l . |22j and 
Voituriez et al. [2| have considered the hydrodynamics of an active polymer solution including 
explicitly the flow of the solvent. In their formulation activity enters via a chemical potential 
proportional to ATP concentration. In our approach this corresponds to the density m of active 



motors. The polarization equation considered in 



3J| contains a term like our ~ Pdjp in Eq. (|5.22|) 



and is obtained there by allowing a coupling between density and splay deformations in the free 
energy of the system. In equilibrium polar fluid this term is ultimately responsible for the instability 
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4C|. 



of a uniformly polarized phase to splay deformation 

A nonequilibrium convective-type term of the form contained on the left hand side of Eq. (|5.22() 
was included in the hydrodynamic equations introduced by Simha and Ramaswamy 3] to describe 
the dynamics of self propelled nematic particles in a solution. In that case the effect of self- 
propulsion was incorporated by assuming that the particles have a mean drift in the direction of 
polarization relative to the solvent, yielding an advection term of the type obtained here. 

In the polarization equation it is apparent that rotational effects from polar crosslinks ('jp) 
increase the bend stiffness, but do not renormalize the splay elastic constant. Nonpolar crosslinks 
(jnp) play a role similar to that of excluded volume corrections in suppressing rotational diffusion. 
This is not surprising as nonpolar crosslinks enhance nematic order in the system. Filament 
bundling described by a renormalizes both the splay and bend stiffness and promotes spatial 
inhomogeneities in the polarization. 



1. Linear Stability 

To examine the stability of the polarized state we choose the y axis along the direction of broken 
symmetry and expand the hydrodynamic fields about their equilibrium values as 



p(r) = po + Sp(r) , 
p(r) =y + 5p ± (r) , 



(5.24) 
(5.25) 



where 5p = x.5p x + 0((5p x ) 2 ). Expanding the fluctuations in Fourier components, the linearized 
equations are given by 



D^k 2 — afhpoky — liwk 



Pk + 



iwk x + ampok x ky 



PoPk 



(5.26) 
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d t pk = - K a (k)k 2 + iw'ky p k + iw"k x - D' a k x k y — , (5-27) 
L J L J p 

where 

w = ^rh(3 Po , (5.28) 
36 

w' = ^rhpp , (5.29) 
3d 

13 

w " = n™PP<> ' ( 5 - 3 °) 
3d 

DZ = j(p-2ampo), (5.31) 

1 / Tfi \ 

D' a = - [D + -^JnpPo - Gampoj , (5.32) 

- 1 /5 77i \ ky 1 / 7 rh \k 2 , 

K a (k) = - y-D + — jppo - mapo J + - y-D + — -f P p - 3map J . (5.33) 

Note that K a (k), with k = k//c, is a generalized stiffness for splay (k y = 0) and bend (k x = 0) 
deformations. Denoting by <fi the angle between the wavevector k and the direction of broken 
symmetry the hydrodynamic modes describing the decay of density and polarization fluctuations 
are given by 

z±(k, 4>) = ikv±(<p) - V^{(f>)k 2 . (5.34) 
The modes are always propagating with speed 



v ±{4>) = {w — w' /2) cos cp ± y (w + w' /2) 2 cos 2 (f> + ww" sin 2 <f> . (5.35) 

The angular dependence of the speed of propagation is shown in Fig. H3 The decay rate is given by 
1 r 



^(4>) - 2 



D a + K a ((j)) — fhapo cos 2 



1 ,{w + w'/2)[D^ - K a ((f)) - rhapo cos 2 <p] + sin 2 4>(wD' a - w"rnap ) ,_ oc , 

±-cos<p .(5.3b) 

^ V (w + w' /2) 2 cos 2 + urn/' sin 

For large values of the bundling rate a the various elastic constants are driven to zero and 2?± (0) < 
0, corresponding to the instability of the uniform polarized state. The condition 2?±(</>) = defines 
the value /?b(0) of the filament density above which the polarized state is unstable. This value 
is largest for <j) = n/2, corresponding to fluctuations with k normal to the direction of mean 
polarization (i.e., pure splay deformations of the local polarization) due to the stiffening of the 
splay elastic constant from polar crosslinks. In contrast, the bend stiffness is not renormalized by 
polar crosslinks, resulting in a lower value of ps at <j) = 0, where polarization deformations are 
pure bend. The angular dependence of /?b(0) is shown in Fig. [7J 



25 



Finally, it is useful to consider explicitly the two limiting cases <f> = and <p = tt/2. For = 
(i.e., k y = k) density and polarization (in this case bend deformations) fluctuations are decoupled. 
Their respective relaxation rates are given by 

r3 5 i 

z + (k,(f>) = z p (k,(f)) = 2iwk — —D — -fhapQ k 2 , (5.37) 

Z-(k, 4>) = z p (k, 4>) = —iw'k — - -D + —jppo — mapo k 2 . (5.38) 

4 L2 6 

The bundling instability is controlled by the growth of density fluctuations and occurs at 

•XT) 

«^ = 0) = Io^- (5 - 39) 

For <j> = ir/2 (i.e., k x = k) the modes are complex conjugate, with 

A' 2 



z±(k, (f> = tt/2) = ±ikVww" - — (13D/2 + rh-/ P po/6 - 9map ) . (5.40) 

8 

Both density and splay fluctuations of the polarization field go unstable at the same density, given 
by 

The zipping effect described by jp tends to stabilize the system. 



D. Summary 

All homogeneous states are rendered unstable by the same mechanism of filament bundling, 
driven by the parameter a. Up to numerical constants and assuming < a, the density above 
which the homogeneous states are unstable can be estimated as pp ~ D/rha. The bundling 
instability line is shown in Fig. |S| One important observation is that the nature of the instability 
changes from diffusive in both the isotropic and the nematic states to oscillatory in the polarized 
state. This suggests that at high filament and motor density the uniform polarized state may be 
replaced by spatially inhomogeneous oscillatory structures such as vortices. 



VI. DISCUSSION 



Several other authors have recently put forward descriptions of the dynamics of active solutions 
and gels of long biofilaments and molecular motors. It is useful to compare our work to others in 
some detail. 
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FIG. 6: The propagation speeds v+(cj>) (blue curve online) and —v^{4>) (red curve online) of the hydrody- 
namic modes in the polarized state as a function of the angle <j). The speed v± is measured in units of lD r 
and we have used m = 1, f3 = 1 and po = pn- 



Kruse et al. 
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22 1 and Voituriez et al. 



3a | have developed a continuum phenomenological de- 



scriptions of the polarized state of active polymer solutions where the hydrodynamic equations are 
written down on the basis of general symmetry considerations. Our work, in contrast, derives such 
equations from a systematic coarse-graining of a more microscopic kinetic equation. The advantage 
of the former method is its generality. The disadvantage is that the resulting continuum equations 
contain many undetermined parameters. Our work yields an estimate for these parameters and an 
understanding for the microscopic mechanisms that control each term in the continuum equations. 
On the other hand, the precise dependence of the parameters on the physical properties of the 
crosslinkers is determined by the specific microscopic model considered, as shown in Appendix lAl 
The two approaches are clearly complementary and both provide insight in the system's dynamics. 



The work described in Refs. 
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2211 and 



3a | explicitly incorporates the dynamics of the solvent, 



which is assumed quiescent in our work (see, however, Ref. [4l|). but consider systems near 
equilibrium by only keeping terms of first order in the chemical potential which controls the rate 
of ATP consumption. For a more precise comparison we refer to Ref. [3^ |. where the equations are 
written in the simpler form appropriate to an active viscous solution, with no viscoelastic effects. In 
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FIG. 7: (color online) The critical density pg where the homogeneous polarized state becomes unstable 
versus angle <p for D = 1, jp = 1, jnp = 1- As in the nematic phase, by increasing a or rh the stable 
region decreases. The critical density is not very sensitive to changes in jp. 



Ref. 38] activity is controlled by the difference A/i in chemical potential of ATP and its hydrolysis 
products, assumed to be constant. This corresponds in our work to the product mug = maRATP 



1 active contributions are proportional to the 
corresponds to our polar rotational rate 



which controls ATP consumption in the system. A 
combinations ttiRatp- The parameter A of Ref. 
7p, describing the "zipping" of filaments due to the action of polar crosslinkers and responsible for 
establishing the homogeneous polarized state. This term is ignored in Ref. [3^ |. where it is assumed 
from the outset that the system exists in a polarized state with pn = 1. The terms proportional 
to the parameter w arising in the polarization equation of Ref. [38j have the same structure as the 
first term on the rig ht hand side of our polarization equation, Eq. (|5.22j) . However, this term is 
obtained in Ref. |3Ja] as derivative of a phenomenological free energy, consisting of the usual Frank 
free energy for a nematic plus a term ~ uipV • p, which is allowed in a polar fluid. As a result, 
this term does not appear explicitly as an active term proportional to A//. On the other hand, 
our analysis of the homogeneous states shows that activity is probably crucial for establishing the 
polarized state as the zipping rate 7p induced by polar crosslinkers is likely to depend on ATP 
consumption rate. Furthermore, the convective terms on the left hand side of Eqs. (|5.21|) and ()5.22|) 
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FIG. 8: Bundling renders homogeneous states unstable for po > ps, where they are replaced by inhomoge- 
neous solutions. The ps line may lie above the Pnp — Pip line or cross through the N and /states, as shown 
in the figure ('Jp/jnp = 1) c% = 50, 7jvp/« = 0.6), depending on the value of Jnp/cx. The instability of 
the / and iV states is diffusive (dashed line), while the instability of the P state is oscillatory (dotted line). 

cannot be obtained as derivatives of a free energy and are therefore not present in Ref. jssl ]. These 
terms are linear in the ATP consumption rate and are important in controlling the oscillatory 
nature of the spatial patterns (e.g., vortices) that are obtained at high filament and motor density 
H- 

Aranson and Tsimring 



3 3l ] have used a generalization of the Maxwell model of binary collisions 
in a gas to describe the dynamics of a solution of polar rods with inelastic interaction representing 
the effect of active crosslinkers. Although their kinetic model, in contrast to ours, allows for instan- 
taneous large changes in the relative angle of two rods upon collision, the continuum equations for 
density and polarization obtained from the model have the same structure as ours. Our parameter 
a corresponds essentially to their parameter B 2 (related to the spatial range of the interaction 
between two rods), while our parameter (5 is proportional to their parameter H, which controls the 
strength of the dependence of the interaction between two rods on their relative orientation (al- 
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though again these authors do not include the convective terms ~ in the density and polarization 
equations). The dependence on motor density or ATP consumption rate does not appear explicitly 
in the continuum equations of Ref. j^] as the strength of the motor-mediated interactions is scaled 
out of the calculation. 

One important difference between our work and that of Ref. [3| is that by incorporating 
excluded volume effects and including the action of both stationary and mobile crosslinkers, we 
can obtain a complete characterization of the homogeneous states of the system. In particular, 
we show that both nematic and polar order are possible in different regions of parameters, and 
evaluate the effect of crosslinks on the isotropic-nematic transition. 

Our work can be extended in several ways. First, we have assumed that the solvent is quiescent 
and only provides the damping on the dynamics of filaments and motors. Relaxing this approxi- 
mation requires considering explicitly the dynamics of a two-component system. In particular the 
dynamics of the solvent must be incorporated when considering the response of the system to an 
externally imposed flow. This will be discussed in a future publication |4lJ. Secondly, an analysis 
of the nonlinear equations for the director and polarization fields in the nematic and polarized 
phases, respectively, reveals the structure of the possible topological defects in each phase and 
their stability. This analysis can be carried out partly analytically and partly numerically and 
can be used to study the range of stability of the spatially inhomogeneous patterns seen in the 
in vitro experiments by considering each pattern as composed of topological defects of the bulk 
system. Finally, for comparison with experimental systems it is crucial to consider the dynamics 
of active solutions in specific geometries, with suitable boundary conditions An important 

application of the dynamics of active filament solutions and gels is that of cell locomotion on a 
substrate. This may be modeled by considering a thin active layer on a substrate, but will require 
incorporating in the model the nonequilibrium polymerization-depolymerization of the filaments, 
mechanical coupling to the substrate and understanding the interplay between them and activity. 



Acknowledgments 

We thank Sriram Ramaswamy for useful discussions. This work was supported by the Royal 
Society (TBL) and by the National Science Foundation, grants DMR-0305407 and DMR-0219292 
(AA and MCM). 



30 



APPENDIX A: MICROSCOPIC MODELS OF MOTOR-FILAMENT KINEMATIC 

In this Appendix we describe some microscopic models of the motor-mediated interaction among 
two filaments. Clearly such models are a great simplification of the contributions to the motor- 
mediated forces, but they allow us to estimate the various phenomenological parameters introduced 
in Section [H] and to justify the approximations used in this paper. We consider two classes of 
models: 1) small motor clusters with an inhomogeneous stepping velocity that vanishes at the 
plus end of the filament, inspired by kinesin motor constructs interacting with microtubules; 2) 
filamentous motor clusters with an antiparallel arrangement of heads inspired by thick myosin 
filaments interacting with thin actin filaments. 



1. Stalling clusters 



The model presented here extends the one discussed in Ref. |2J] to include the possibility of 
motor-induced filament rotation. 

We consider a pair of filaments (denoted as filaments 1 and 2) cross-linked by a motor cluster. 
Due to the action of the motors, filaments 1 and 2 acquire center-of-mass velocities vi and V2 and 
rotational velocities u\ and u>2 about the center of mass. Our goal is to evaluate these velocities in 
terms of the rates at which the motor cluster steps along the filament and rotates relative to it, and 
of the filaments' orientation. Both filaments and motors move through a solution. We assume that 
the filament dynamics is overdamped and the friction of motors is very small compared to that of 
filaments. The coupling of the filaments to the fluid is via a local friction (Rouse model). This 
is a reasonable approximation for a quiescent solution without externally imposed flow nor flow 
generated by the motor activity. Under these conditions, hydrodynamic coupling yields logarithmic 



corrections to the friction, which are small for long thin rods |9|, llC 



)dyna 

Id, Lo 



271 ] . Momentum conservation 



requires that in the absence of external forces and torques, the total force (torque) acting on 
filaments centered at a given position be balanced by the frictional force (torque) experienced 
by the filament while moving through the fluid. The third law requires that any force or torque 
generated by an active crosslink on one of the filaments of the pair is balanced by an equal and 
opposite force or torque acting on the other filament. This yields 

Cy(ttl)«lj = -Qij{u2)V2j , (Al) 
Cr^l = ~( r U2 , (A2) 

where dj(u) = (\\UiUj+(±(5ij—UiUj) is the friction tensor of the rod, with (y and (± the longitudinal 
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FIG. 9: Two filaments of orientation u\ and U2 connected by an active torsional spring. Here 9i = 6(si) , 
for i = 1, 2 is the torsional angle at the point of attachment of the motor cluster. Note that we have chosen 
a convention such that 6 = Q\ — #2- 



and transverse friction coefficients, respectively, and £ r is the rotational friction. Equation (|Aljl 
shows that the anisotropy of the friction tensor allows for a net translation (vi + V2 7^ 0) of the 
filament pair induced by motors. 

The mobile crosslink is a small aggregate of molecular motors that exerts forces and torques 
on the filaments by converting chemical energy from the hydrolysis of ATP into mechanical work. 
While walking along the filaments, the motor clusters can also apply aligning torques on the 
filaments, if there is a preferred angle between the heads of the motor cluster. To capture this, we 
consider the cluster to be a nonlinear torsional spring of size l m ~ b « I. A similar description 
would also be appropriate for any polar cross-linking protein. However, a motor cluster which 
aligns the filaments by active contractions has an ATP-dependent spring constant. A schematic is 
shown in Fig. El 

It is convenient to think of the motor cluster as composed of two heads, with the i-th head 
(i = 1,2) attached to the i-th filament at position r* = + MjSj. Motor heads are assumed to 
step towards the polar end of filaments at a known speed, u(s), which depends on the point of 
attachment (see Fig. 110(1 . The motor-induced torques occur along the axis of the motor cluster, 
assumed to be directed perpendicular to the plane containing the filaments, and are capable of 
generating equal and opposite torques on the two filaments. The torsional angles 8(s) obey the 
following equations Cr^(si) = — Cr@{s2) = — « sin [Q(s\)— 9{s2)\ — —k[0(si)—9(s2)], k is the torsional 
spring constant. In general the torsional spring constant will also depend on the position of the 
motor cluster along the filament, i.e., k = k(s). The resulting inhomogeneities in the rotational 
rate 7 does not yield qualitatively new terms in the hydrodynamic equations and will be neglected 



here (see also 



37 1). The dynamics of the i-th motor head is described by a translational velocity 
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FIG. 10: (a) Two filaments connected by an active cross-link and the geometry of the overlap. The filaments' 
centers are separated by £ = sitti — S2U2- (b) The profile of the motor stepping rate. 

v™ at the point of attachment and a rotational velocity u>™, given by 

v™ = if = Vi + u(s,)tii + x ttj , (A3) 
+ (A4) 

C r |t*l X U 2 | 

with u(sj) = Sj . Finally, the two motors within a cluster are rigidly attached to each other. This 
requires 



'2 > 



(A5) 
(A6) 



Since the motor cluster has size Z m ~ b « I, we can neglect the length of the motor compared 
to that of the filament and assume that the attachment points satisfy r * = r 2 , or £ = r 2 — ri = 
siUi — s 2 w 2 . Equations (|A1[) and (|A2|) . together with the expressions (|A3|) and (|A4j) for the 
velocities of the motor heads, and the conditions ()A5|) and (|A6|) that the two motor heads are 
rigidly connected, then provide a closed set of equations that can be solved to obtain the filaments' 
translational and rotational velocities in terms of their relative orientations and of the motors' 
stepping and torsion rates. It is convenient to introduce relative and net translational and rotational 
velocities of the filament pair as 



V 



vi - v 2 , 
Vi + v 2 



(A7) 
(A8) 



and 



to 



U>! - CJ 2 . 



(A9) 
(A10) 
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For the translational velocities we obtain 



v = u(s 2 )u 2 - u(s\)ui + 7(s 2 - si)(«i + u 2 ) + 7(1 + Ul ■ U2)(siUl - S 2 U 2 ) 
= ^|[n(s 2 ) - u(si)} - 7(si - s 2 )(l - ui ■ u 2 )}{u 2 + ui) 

+^{N s 2) + w(si)] - 7(si + s 2 )(l + «l • U 2 )}(u 2 - Ul) , 



(All) 



with 



V = 4(u 2 + ui) + B(u 2 - ui) , 



,4 
D 



a / 1 — ui • u 2 

4 V 1 — (TUi • u 2 
cr / 1 + Ui ■ U 2 



4 V 1 + aii\ • u 2 



u(s 2 ) + u(si) - j(si + s 2 )(l + Ui • u 2 
u(s 2 ) - u(si) + 7(s 2 - si)(l - ui • u 2 



(A12) 

(A13) 
(A14) 



where we have defined 7 = k/C t and a = 7-^ — ^4. For long thin rods C± = 2Cn = 2C and cr = 1/3. 

(C± + C||) " 
There is no net rotational velocity of the pair (fi = 0). The relative rotational velocity is given 

by 



Ui X u 2 . 

w = 27 —± 7—7 . (A15 

|ui X w 2 | 

The fact that V ^ indicates that motor activity can induce a net motion of the pair relative to 
the solution. This is a consequence of hydrodynamics at low Reynolds numbers, which gives an 
anisotropic of friction tensor for long thin rods. As a result V vanishes when (j_ = £||- Also V 
vanishes identically for u 2 = ±ui, so that V = in one dimension. 

We can compare the expression for the filament velocities obtained via the microscopic model de- 
scribed in this section to the general expression introduced on the basis of symmetry considerations 
in Eqs. (|2.1UB2.11() by expanding the stepping rate as u(s) « Uq — su', where v! = —du{s)/ds > 0. 
Substituting the expressions ^[u(si) — u(s 2 )] ~ -y ( Sl — s 2 ) and ^[u(s 2 ) + u(si) ~ uq + \{s 2 + si) 
into Eq. (|A11|) we obtain a general expression for the relative velocity given by 

v = a + u + (£ • u+) + a_u_(£ • u_) + (3(u 2 - ui) , (A16) 

where u + = (u 2 + ui)/|-u 2 + ui| and u_ = (u 2 — ui)/|u 2 — ui| and 

a + = -7(1 - ui • u 2 ) + u' , (A17) 
a_ = 7(1 + ui -u 2 ) +u' , (A18) 
P = u . (A19) 
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If u' 3> 7, then a+ = a_ = a, leading to the simpler expression for the relative velocity (see Eqs. 
(|2.1()I2.11|) ) which we use for the whole of this article . 

By comparing Eqs. (|A16|) . (j A 13(1 and ()A15|) to the general expressions given in Sec. [H]we obtain 
the following estimates 

ao^lA, (A20) 
as 

(3 ~u , (A21) 
IP ~ K/C r . (A22) 

Note that the specific microscopic model used here gives &\ = 0,/?i = and "fNP = 0. This is the 
result of having considered the kinematics of a single pair of filaments coupled by one motor cluster. 
An additional dependence of the effective coupling constants is introduced by the dependence of 
the motor binding probability on the relative orientation of the filaments. 



2. Contractile motor filaments 



Here we consider another microscopic model relevant to large contractile filaments of myosin II 
(thick mini-filaments) interacting with filamentous actin (thin filaments) . Both the thick contractile 
motor filament and the thin filaments undergo overdamped motion in a quiescent fluid. 

We consider two thin (e.g., actin) polar filaments of length I with centers of mass at ri and 
r2 and orientations ii\ and «2, respectively. The active crosslink is a (thick) filament of motor 
proteins of length l m < I. The motor heads within the motor filament are antialigned, with the 
motor heads at the two ends of the motor filament pointing in opposite directions. Its orientation 
is described by a unit vector u m oriented along its axis. We choose the direction of u m to be 
from thin filament 1 to 2, as indicated in Fig. 1111 The motor filament exerts torques on the actin 
filaments by acting as a torsional spring of strength k. As a result, the actin filaments align with 
the motor filament, as shown in Fig. 1111 Once the thin filaments are parallel to the motor filament, 
the heads on either side of the motor filament pull the thin filaments together until they overlap 
over a length l m < I. In this configuration both thin filaments are linked by both heads at the two 
ends of the contractile motor filament. The effect of these two heads balance and the thin filaments 
remain stationary relative to each other. 

To describe the dynamics, we denote by v m and u m the center of mass and angular velocity 
of the thick motor filament, respectively. The friction tensor of the motor filament is given by 
Cij = CmW^miUmj + (m±((>ij — UmiU m j). Since the thick motor filament is shorter than the two 



35 




FIG. 11: Two thin filaments of length I connected by a contractile thick motor filament of length l m . Due 
to the torsional springs, the thick filament aligns the actin filaments in an antiparallcl configuration. The 
'stepping' of motor heads towards the plus ends of the thin filaments bring the (almost) antiparallel filaments 
together if their centres are more than I — l rn apart. When their centers of mass are separated by I — l m , the 
two actin filaments are stationary due to the opposing effects of the two motor heads at the opposite ends 
of the motor filament. 

thin filaments it crosslinks, we expect the motor filament translational friction coefficients £ m y and 
Cm± an d rotational friction Cmr, to be smaller than the corresponding parameters for the actin 
filaments, i.e., Cm||>C||) Cm± < Cu an d Cmr < Cr- The separation of centers of mass of the thin 
filaments is £ = r 2 — ri = u±si — u 2 s 2 + u m l m . Force and torque balance require 

Qij{ui)vij + Cij(u 2 )v 2 j + (%j(Um)Vmj = , (A23) 
Cr<^l = Gl , 
Cr&2 = G2 , 

Cmr^m = -Gl - G 2 , (A24) 

where 

Gi = K(u m x iii) , 

G 2 = -K(um x u 2 ) . (A25) 

The position of the center of mass of the thick motor filament is r m = t\ + s\U\ + u m l m /2 = 
r 2 + 52^2 — u m l m /2 and its velocity is given by 

v™- = r m = v» + u(si)ui + SiUJi xui + (-iy~ ll -^u m x u m , % = 1, 2 . (A26) 
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The set of Eqs. (|A"26jl . (|A"23l) and (|A"24|) can be solved for the mean V = (vi + v 2 )/2 , Q = 
(uii + u 2 ) /2 and relative v = vi — v 2 , u> = u\ — u> 2 translational and angular velocities of the two 
thin filaments and the velocity v m of the thick motor filament. The general solution is complicated 
and not terribly transparent. 

To simplify the solution, we assume that l m <C I so that Q mr <C ( r and £ m *C C- I n this case 
the orientation of the thick motor filament relaxes much faster than that of the two thin filaments 
and is therefore slaved to the thin filament orientations, so that u m = (u 2 - ui)/\u 2 - ui\. The 
expressions for the relative and mean velocities are then 



r \ * 1 \» 7 r(l + W2 • Ml) . w . 

v = u(s 2 )u 2 - u(si)ui+ -< —z r-r- \ u 2 ~ ui) (s 2 + s ly 

2 I \u 2 — Ui\ 



(1 - Ml • U 2 ) 
\u 2 - Ul\ 



(si - s 2 ) (u 2 + tii) } 



(A27) 



V = A (u 2 + ui) + B (u 2 - ui) 



(A28) 



where 



.4 
B 



1 — Ui-U 2 

1 — au\ ■ u 2 
1 + Ui ■ u 2 
1 + au\ ■ u 2 



(1 + u 2 ■ Ui) 
u(s 2 ) + u(si) + 7 — 7—r- [S2 + si, 



u(s 2 ) - u(si) + 7 



\U 2 - Ml| 

[l-ui- u 2 ) 
\u 2 - tii| 



{s\ - s 2 ) 



(A29) 
(A30) 



The center of mass velocity of the thick motor filament is given by 



v m = V + w/2 , 



(A31) 



where 



w = u(s 2 )u 2 + n(si)t*i- -< —rz r-r- («2 - ™i («i - $2) 

2 I ti2 — tti 



+ (1 m-ua) (si + ^ + (A32) 

\U 2 - till 



with a = — There is no net rotational velocity of the pair (fi = 0) and the relative 

(Cl + C||) 
rotational velocity is given by 



c = 2 7 " 2 X . (A33) 
t*i x «2 



where 7 = 
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The speed of the motor clusters depends on the thin filaments' relative position and orientation: 
it is maximal when the filaments are oppositely oriented and goes to zero when the filament overlaps 
a length l m (see Fig. Illj) . A simple expression which satisfies these conditions is 

1 I, 



(A34) 



where uq is the (constant) single-motor 'stepping rate' (step-size/cycle-time). 

Assuming uq 3> 7 and substituting Eq. QA34j) into Eqs. (|A27[) . (jA28j) and (|A33|) . we obtain an 
expression for relative velocity of the two thin filaments connected by a thick filament, 

Z 1(1- lm) 



V ~ Uq(1 — U\ ■ U2) 



(U 2 ~ Ml) 



(A35) 



I 2 I 

where £ = siiii — S2U2 + l m u m . In obtaining Eq. ()A35|) we used approximated u m ~ 1/2(142 — 111), 



which holds when the filaments are antiparallel. By comparing Eqs. (|A35|) and ()A33|) to the general 
expressions given in Sec. [H]we obtain the following estimates 



a 



A) = -Pi -j-) u ° ' 

IP = 7 - «/Cr • 



(A36) 
(A37) 
(A38) 



The crucial difference between the effect of the stalling crosslinker considered in Sec. IA II and the 
contractile minifilament considered here is of course that in the present case the rates do (/3o) an d 
Si (/?i) have opposite signs. This will have important effects on the system's rheology. 



APPENDIX B: MOMENT EXPANSION 



To define the exact moment expansion of the filament concentration c(r, u, t) we introduce a set 
of irreducible tensors i which are equivalent to the spherical harmonics, but are expressed 

in Cartesian coordinates. The components of i are homogeneous polynomials of degree 

m in the components of the unit vector u, with the properties that they are fully symmetric in 
the subscripts i\, 12, i m , and that no nonvanishing tensor of lower order can be formed by 
contraction. 

Here we denote by d the dimensionality and write the general expression of the T m for d = 2,3. 
Each tensor is orthogonal to all the other ones and normalized according to a product defined by 

(a,b) = J ^ <.(«)&(«), (Bl) 
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where Qd is the solid angle in d dimensions, with = 27r and $7 3 = Am. The first few irreducible 
tensors are 



T° = l 



(B2) 



T} = u t , 



(B3) 



Uiilj Sij , 



(B4) 



„3 _ „ „ „ 

J jjfc — u i u j u k — ^ — 2 



(B5) 



ji4 



1 



ijfcZ ~~ u i u j u k u l 



d(d + 2) 



(B6) 



Ujfcis = UiUjU k UlU s - 



{d + 2){d + 4) 



UiA jkls + UjA 



j^ikls 



Uk\jls + Ul^ijks + U s Aij k i 



(B7) 



where repeated indices are summed over and the tensor Ay-jy is given by 



Aijki = 5ij5ki + SikSji + SuSjk 



(B8) 



it also satisfies A^kk = + 2). 

The filament concentration has an exact expansion on the basis of these irreducible tensors, 
given by 



c(r>,t)=^< 42 ... m (r,t)^ 2 .., m (^) 



m=0 

where the m-th order moment a™„ ,■_ is a tensor determined by 



1112... %m 



m 



1 



Q hi2-im 



The first three moments are given by 



(B9) 



(u)c(r,u,t) . (BIO) 



du „ 1 . 

— c(r,u,t) = — p(r,t) , 

a i = d [ ?r u% c(r,u,t) = —p(r,t)pi(r,t) 



a?- = d ( d + 2 ) y ^ (niUj - ^)c(r,u,i) = d '^" 2 \ (r, t)S i:j (r, t) . 



1.) 



(Bll) 
(B12) 
(B13) 



39 



Here p, p and Sij are the density, polarization, and nematic order parameter of the rods, respec- 
tively. Retaining only moments up to the third one, the filament concentration can be written 

as 



1 



c(r, u, t) « t p(r, t) 1 + du- p(r, t) + 



d(d + 2) 



Qij(u)Sij(r,t) 



(B14) 



with Qij = UiUj - \S 



i.r 



APPENDIX C: DERIVATION OF COARSE-GRAINED CURRENTS 



In this Appendix we outline the derivation of the currents and source terms entering the equa- 
tions for the density, polarization and alignment tensor and give their general form. For simplicity 
we restrict ourselves to the case of long thin rods, where D\\ = 2D± = D and a = 1/3. 

It is convenient to separate the translational and rotational currents defined in Eqs. (|3.8B3.1U|) 
and (|3. 1 1113.12]) in diffusive, excluded volume and active contributions, 

J% = J{ + </j + J$ , 



Jij - 
Jijk 



tP 4- T ex + / 1 

l J l 3 l 3 ' 

- T D 4- f cx 4- T A 

~ J ijk * J ijk * J ijk ' 



and 



T> TfD I OCX , J~>A 

rLi — -ft,- ■+■ rt,- -f- it, 



R 



R 



i i 

A 
'J 



(CI) 

(C2) 
(C3) 



where each contribution arises from the corresponding term in Eqs. (|2.2B2.3|) . 

The diffusive contributions are evaluated by inserting the truncated moment expansion for 
the filament concentration in the corresponding contributions to the translational and rotational 
diffusion currents in Eqs. ()2.2II2.3|) and performing the angular average, with the result 



J 



D 
ij 

D 



a 



D 



D 



a 



D 



_D(?> 

5 



ijk 



, Su v • t + -dm 

4 V J 2 3 

[Sjkdip + S ik djp - 5ijd k p) - — 



(C4) 
(C5) 



SikdiQji + SjkdiQu + -dkQij - SijdiQki , (C6) 



and 



R 



D 



D r Ti , 



Rn — ^D r Qij . 



(C7) 
(C8) 
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To evaluate the excluded volume contributions we expand the concentration in Eq. ()2.4|) near 
its value at ri as in Eq. (j3.2j) . truncate the moment expansion of the concentration to third order 
and perform the angular integrations. Retaining terms up to first order in the gradients of the 
fields in the currents and up to second order in the source terms, we obtain 



J? 



Dv dj 



-p(Qij + ^ijpj + gQikQjk + Y^^ijQkiQki 



+ -DvopdjQij , (C9) 



TCX 

J ij 



-Dv 
2 



-(Ajjjy + A5u5jk)T k dip - -(Aijin + 6S in 5ji)T k d n QM 



(CIO) 



where vq = -, and 



J ijk = Y2 DV °\ ~ g( SikS i l + SilS i k ~ Si i Skl "> dl P 

7$kiQij + 5ikQji + SjiQik + SuQjk + SjkQu — Z&ijQki 



75kiQij + SikQji + SjiQik + SuQjk + bjkQu — 2$ijQki dip 

QijQkl QikQjl QilQjk ^ijQkrQlr 
+ (5uQj r + 5jlQi r — Sij 

Qlr)Qkr ~i~ ($ikQjr ~t~ ^jkQir $ijQkr)Qlr 
+ -^(3ik5jl + Sildjk — 198ij5kl)QrsQrs + ^^klQirQjr j , 



D r v 



6 
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+ ( TfrdkdjQij + TjdidkQjk — TidjdkQjk 



(Cll) 



(C12) 



rpex 



-D r v pQ 



i.i 



D r v 



p[didjp - -SijV 2 p) + (5uQjk + 5jiQ ik - SijQki)dkdip (C13) 



To evaluate the active contributions, we insert the gradient expansion of the concentration and 
motor density given in Eqs. (|3.2|) and 1)3. 3|) in the Eq. 1)2.8)1 for the motor current density and in 
Eqs. (|2.5|) and ()2.6)1 for the filament currents. The integrals over the lengths s\ and s 2 of the 
filaments can then be evaluated explicitly. All terms containing odd powers of components of the 
filament center of mass separation £ vanish when averaged over the rods' length. To evaluate the 
angular integrals in the filamment current densities we also expand the translational and rates a{9) 



and (3(9), as well as the excluded volume x u 2 \ = a/1 — («i • U2), to first order in the cosine of 
the angle between the two filaments, u\ ■ U2- With this approximation, the motor- induced linear 
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and angular velocities vi = V + v/2 and u\ are written as 



y(U 2 - Ml) + —i 



+ 



^(£ 2 - u x ) + ^] (™i • 62) + £>((«! • «2) 2 ) 



(C14) 



V = ~flo(ui + U2) + ^ao{u\Si + u 2 s 2 ) 
Po(l - cr)](ui • M2)(«l + w 2 ) 
+^[5i(wisi + m 2 s 2 ) - qq(1 - cr)(wis 2 + m 2 si)](wi • tt 2 ) + 0((u\ • 'U 2 ) 2 ) , (C15) 



u>i = 2[7p + 7tvp(wi • w 2 )]'Ui x ti 2 + 0((u\ ■ u 2 ) 



(C16) 



and |wi X tt 2 | ~ 1. As indicated in the main text, contributions of higher order in Ui -it 2 only change 
the values of the numerical coefficients of the various terms in the expressions for the currents given 
below, but do not contribute any qualitatively new terms. 

Finally, we insert the moment expansion of the filament concentration c(r, &,t), truncate it to 



the first three moments, as given in Eq. (|B14|) . and evaluate the active contributions to the various 
current densities defined in Eqs. (|3.8W3.1Uj) and (|3.11B3.12() . The calculation of the angular integrals 
is quite lengthy and has been carried out with Maple. 
The motor current density is given by 



jm 



48 



Tjdidjm 



+ C(V 3 ) . 



(C17) 



The active contribution to the current density is naturally written as the sum of two parts 



J^(T,t) = pV i + d j af j , 



(C18) 



where 



pV 4 



+ - 



1 



Qij + 2^ P ) ~ \ a ° ( ® ik + ^ikP) 



djm , 



and the active contribution to the stress tensor, cr^, is given by 



1 



1 



(C19) 



(C20) 



with ao = ao/48 and a\ = di/48. The drift Vi vanishes in a passive system and arises entirely from 
the contribution to the active current from the net velocity V of the pair. It is in fact proportional 
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to a = — C||)/(C± + Cll) an d vanishes for isotropic objects. The term proportional to ceo m the 
stress tensor describes the built-up of density inhomogeneities via filament bundling and has an 
effect opposite to that of conventional diffusion. As shown below, this is the main term responsible 
for driving the instability of homogeneous states. 

The active contributions to the translational and rotational polarization currents are given by 



where 



and 



J; 



6 3 



m 
~6 



(Pi ~ A)/6) 



T-Tj + \k-^ 



-m^-p(Qij + + y(A + 0o/3)(Qik + ^kkpj (Qjk + \^jkP 

+^-m(T i d j + Tjdi + S i:i T • V)p + ^rh P T tj + ^rhT^p 

+^-rhT j d k Qik + 5ijTid k Q k i - T k d k Qij + Tid k Q kj - T k djQ ik 

+-^-rnQj k T ki H — -^-fh -5ijQ k iT k i — QijV ■ T + Qi k T k j + Qj k T ki 

+-(«! + lla /6)p(T i d j + Tjdi + SijT ■ V)m 



+- 



2ai 
9~L 



Q jk (T k di + Tid k ) + 5 l] Q kl T l d k 



m 



+-(2«! - Q /2) Qik(Tjd k + T k dj) + Q l3 T k d k 



m 



(C21) 



= da , j + d j T i + 5 ii V-T . 



(C22) 



--fpfhpTi + (27p - j N p)rhTjQ 



TP rl 

24 1 4 
1 



p(3<%V 2 - 2d i d j )(rhT j 



3 L 



Q^(^fcV 2 + 2d j d k )(mT k ) + Qjkdji^diimn) - 5d k (rhTi) 
+^m /0 (<%V 2 + 23 i i )2j - ^mQ^^V 2 + 2djd k )T k } 
~^{\ [ T ^ 6ikV2 ~ d i d k)(rnQjk) 

+2T k d k d J (mQ ij ) - -TidjdkirhQjk)] + -mT^dfi, - ^(%V 2 )p 
+^m[r :? (^ fc V 2 + 2did k )Q jk + 2T k d k djQij - 2Tidjd k Q jk ^ } . 



(C23) 



Finally, the translational and rotational contributions to the alignment tensor current are 



J, 



m 
~Y2 



\ (y A) + Pi) P^ikTj + S jk Ti - 5 i3 T k ) + I (^fl, - 2fc 



1 /5 



+ g (^3^0 _ A ) (TjQik + TiQj k - dijTiQki) 
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i n 



jA) — 2(3i^}Ti(5ikQji + SjkQu — SijQki)} 



+ 



1 19 

-{SikSjn + Sjk&in — 5ij5kn)p + -^$knQij 



u 



3 V3 

+ -(5jnQik + 5ikQnj + hnQjk + SjkQin ~ 2SijQ kn ) 9l(Ql n + -5l n p) 

11 rl 1 



+- 



4 (<5ifc<5j7 + - 5ij5ki)p + -SkiQij 



+^jiQik + hkQji + ^Qjfc + SjkQu - 2$ijQki) di(fhp) 

-^(Sindjl + 5il5j n — 5ij5l n )p + -Sl n Qij 



+7^(duQjn + SjnQu + SjiQi n + 5 in Qji - 25ijQi n ) di(mQ kn )} 
+^ai^-m(S in Tj + 5 jn Ti - 5ijT n )di(5 k iT n + 5 kn Ti + Si n T k ) 
+ 1 -(A mjkln -35 lJ A klnv )T v d l (rhT n )} , 



(C24) 



\jkl — bijSkl + <>i k 5ji + 5u5j k . 



Aijklnp — &ij A k i n p + 5i k Aji n p + 5i[A k j n p -\- 5i n A k ijp + 5ipA k i n j 



(C25) 
(C26) 



-2~f P m(TiTj - ^SijT 2 ) - jNpmpQij 



~ ill?^^ + T ^ 2Tl) + ™ (Ti ^' + T A) d kT k ~ h ijr h(T k V 2 T k + 2T k d k diT{) 
+^ [TiV^mTj) + TjV^mTi) - (T^ + T^d^fhTj,) - T k didj{wT k ) 
-i^(r fc V 2 (mT fc ) + 2T k d k di(fhTi)) + 2T fc 3 fe ($(mTj) + ^(m^))] } 



7^P fm 
" 48 



1 2 

«%<9jP - 2 Si i^ 2 P + g^^i + ZdidkQjk + 2djd k Q ik - 25 ij d k d l Q k i) 



m 

+ 3 
2fh 



Qikdjd k p + Qj k did k p - QijV 2 p - SijQ k id k dip 



Qik^ 2 Qjk + QjkS 2 Qik + 2Qu(djd k Q k i + d k diQj k ) 



+2Q j i(d i d k Q k i + d k diQ ik ) - 4Qijd k diQ ki - 5ijQ kl (V 2 Q k i + 4did r Q 



Ax / 



+ 3" 

1 

+ 3 



2V (rhQij) + did k {mQj k ) + djd k {rhQ ik ) - 5ijd k di(mQ k 



ki, 



QikV 2 (mQjk) + QjkV 2 (rhQ ik ) - 2Qu[d j d k {rhQ k i) - d k di{fhQ jk ) 
-2Qji(did k (mQ k i) - did k (mQ ik )^j + 2Q kl (did k (rhQ + djd k (mQu) 
+3Q k id k di(mQij) - Qijd k di(mQ k i) - QkAdjirnQki) 
-5ijQ k i(^V 2 (mQ kl ) + 2d l d r {rhQ kr fj j . 



(C27) 
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FIG. 12: (color online) The density mode z p as a function of k for different values of the parameter a for 
fh = 1, po = 1, ct c — 0.57 and a c ~ 1.23. For a c < a < a c the isotropic state is always stable, and for 
a > a c the isotropic state is unstable for long wavelength. 

The general nonlinear equations are fairly complicated, but the various terms have simple 
physical interpretations, as will become apparent below. The terms proportional to ao and a\ tend 
to bundle filaments together, therefore enhancing density fluctuations. The terms proportional 
to f3o and (3\ tend to align the filaments in the direction of the polarization and thus suppress 
polarization fluctuations. The jp and jnp terms rotate and align filament and play a crucial role 
in controlling the possible homogeneous states of the system. 



APPENDIX D: ROLE OF HIGHER ORDER GRADIENTS 



Here we discuss the role of terms of order higher than second in the gradients of the hydrody- 
namic fields in controlling the bundling instability. For simplicity we only consider the instability 
of the isotropic state. In this case the only hydrodynamic variable is the filament density. Including 
terms of order up to k 4 , the dynamics of the Fourier components of density fluctuations defined in 
Eq. 1)5. 2|) is governed by 



OtPk 



-k 1 



-D(l + v po) - ampo 



Pk + 



1 

96 



13 



-Dvqpq 



19 



-ampo 



(Dl) 
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Their decay is controlled by a single diffusive mode, given by 

Zp = -C 2 k 2 + C A k 4 , (D2) 

where 



C 2 = mp (a c - a) , (D3) 
19 
480 5 



19 

C4 = -T^ m Po( a 'c ~ a ) > (P4) 



and 



3Dvq 1 + v p 

a c = -r- — — . (D5) 



4m vqPq 
35Dvq 
72m 

At the low filament densities where the isotropic phase exists the value a c where the coefficient 



4 = ^. (D6) 



C 2 changes sign grows rapidly with filament density, while at the value a' c where the coefficient C4 
changes sign is independent of pq. We therefore expect a c > a' c in the region of interest. We can 
then identify three regions: 

• For a < a' c both C 2 and C4 are positive. Long wavelength density fluctuations always decay 
and the isotropic state is stable. The growth rate defined in Eq. (jD2|) becomes positive for 



k > ko, with ko = y C2/C4, but this short scale instability is outside the range of validity of 
the present work. We expect that it will be suppressed by terms of even higher order in the 
gradients. 

For a' c < a < a c we have C 2 > and C4 < and the isotropic state is always stable. 

For a > a c the eigenvalue z p (k) controlling the dynamics of density fluctuations becomes 
positive for k < ko- In this regime long wavelength density fluctuations grow in time and the 
isotropic state is unstable. The isotropic state is stabilized again at short scales, k > ko = 



\J C2/C4. 

The location of the instability in the (a, po) is not affected by terms beyond quadratic in the 
gradients. These terms do, however, introduce a length scale corresponding to the wavevector 
ko = \J C2/C4 ~ \/ (a — a c ) I [a — a' c ) beyond which the isotropic state is stabilized by short scale 
effects as seen in Fig. 1121 The wavevector of the fastest growing mode in the unstable region is 
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k m = \J C2/2C4 ~ \J a c j (a c — a' c ) ~ e 1 / 2 , that vanishes with the distance e = (a — a c )/a c from 
the instability. 
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